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■'  Summary  of  Principal  Recent  Results 

To  date,  the  successes  of  our  effort  for  the  HEDM  program,  both  methodological  and  in  terms 
of  metastable  molecule  chemistry  have  been  many.  Considering  the  latter  first,  recognizing  the 
prospects  for  highly  energetic  metastable  species  based  upon  nitrogen,  we  initiated  our  study 
of  a  variety  of  such  potential  molecules.  These  included  N3H3  [1]*,  which  is  isoelectronic 
with  the  well-known  species  ozone,  cyclopropane  and  propene;  tetrahedral  N4  [6]  which  is 
isovalent  with  P4  or  AS4,  and  is  forbidden  by  symmetry  from  decaying  to  ground  state  N2 
molecules;  octahedral  Ng  [6],  isoelectronic  with  cubane  and  similarly  forbidden  from  direct 
decomposition  to  N2  in  its  ground  state;  and  pentazole,  HN5  [5].  Each  of  these  was  found 
to  contribute  a  metastable  minimum  on  the  global  potential  energy  surface.  The  minima  were 
fully  characterized  by  calculation  of  the  Hessian  (second  derivative  matrix)  frequently  using  our 
recently  developed  analytical  gradient  CC/MBPT  methods,  which  also  provided  the  IR  spectra 
of  these  unknown  molecules.  At  this  level  of  approximation  the  vibrational  frequencies  are 
normally  within  5%  of  experiment  and  the  relative  intensities  typically  to  within  about  20%  [6]. 
I  might  add,  the  character  of  a  critical  point  on  an  energy  surface  can  change  depending  upon 
whether  electron  correlation  is  included,  so  this  is  a  critical  test  for  every  potential  metastable 
molecule.  Furthermore,  some  regard  must  be  paid  to  low-lying  excited  states  that  might  offer 
a  decomposition  path  to  ensure  that  the  molecule  is  stable  to  unimolecular  decomposition.  In 
the  absence  of  matrix  isolation,  bimolecular  mechanisms  can  also  offer  ways  to  prohibit  the 
existence  of  any  such  proposed  molecule. 

Although  sometimes  postulated,  there  is  no  experimental  evidence  for  the  existence  of  any  of 
these  molecules.  This  tends  to  be  a  recurring  theme  in  the  investigation  of  metastable  molecules. 
Yet,  the  accuracy  of.  modem  ab  initic  correlated  predictions  is  sufficiently  reliable  that  barring 
rapid  bimolecular  kinetics  and  certain  possible  excited  state  decomposition  paths,  the  proposed 
metastable  molecule  should  be  capable  of  synthesis  at  least  in  matrix  isolation  and  can  be 
identified  from  our  predicted  frequencies. 


*References  are  to  HEDM  publication  list,  pgs.  15,  16. 
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A.  Metastable  Mdlecules 

1.  The  nitrogen  aromatic  analogue  of  benzene,  hydrogen  pentazole, 
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Figure  1  Hydrogen  Pentazole  Molecule 


has  never  been  observed  experimentally.  Consequently,  it  is  likely  to  be  metastable.  Phenyl 
derivatives  of  pentazole  have  been  synthesized.  We  performed  the  initial  ab  initio  study 
of  HN5[5]  predicting  its  vibrational  spectra  for  future  identification,  and  addressing  its 
decomposition  pathways  and  associated  activation  barriers.  HN5  was  found  to  offer  an 
Isp=346  sec  as  a  monopropellant,  and  the  barrier  to  decomposition  was  found  to  be  19.8 
kcal/mol.  This  study  employed  the  new  MBPT(2)  analytical  gradient  methods  previously 
developed  under  this  grant. 

2.  The  fact  that  HN5  can  be  stabilized  b.y  replacing  H  by  a  phenyl  group  suggests  the  feasibility 
of  using  one  N5  to  stabilize  another,  to  form  the  bipentazole  molecule  [5]. 


Figure  2  Bipentazole  Molecule 

We  believe  that  this  fascinating  molecule  is  an  outstanding  candidate  for  synthesis.  Various 
other  pentazole  derivatives  have  been  considered,  with  amine  (-NH2)  substituted  pentazoles 
(H2N-N5)  showing  greater  Igp  values  than  nitro  O2N-N5  substituted  ones,  but  electron 
donating  substituents  like  NH2  and  Li  also  lower  the  activation  barrier  to  decomposition. 
Electron  withdrawing  groups  stabilize  the  pentazole  ring. 
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3.  Since  Be(Bilt)2  has  one  of  the  highest  heats  of  formation  known,  we  consider  another 
experimentally  unknown  molecule,  LiBe(BH4)3  where  a  third  energetic  BH4  group  could  be 
added  to  the  basic  structure^  [24], 
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Figure  3  Be(BH4)3  Anion 


We  find  that  this  molecule  has  a  substantial  barrier  to  losing  its  third  BH4  unit.  Its  Isp=379 
sec.  Like  bipentazole,  it  deserves  synthesis  and  experimental  characterization. 

4.  Another  group  of  fascinating,  metastable  molecules  are  those  composed  of  HBNCH  units 
[3,7,18].  Our  original  interest  in  this  unit  arose  from  its  isovalent  equivalence  to  the  organic 
analog  cyclopropenyl  cation. 


H 


N 


Figure  4  Cyclopropenyl  Molecule 


We  performed  an  exhaustive  study  of  the  HCNBH  potential  energy  surface  [7].  When  two 
such  units  are  combined,  a  cyclic  molecule  of  the  form 
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Figure  5  Cylic  Dimers  of  HC-N-(BH)-(CH) 

could  be  envisioned,  but  rather  than  being  “benzene-like,”  the  molecule  assumes  a  bent 
form  analogous  to  cyclohexane.  We  present  its  energetics  and  vibrational  spectra  [3,  18]. 
Perhaps  even  more  intriguing  is  the  prospect  for  polymers  of  the  HCNBH  unit  [3,  18]! 


Figure  6  Polymers  of  N-(BH)-(CH) 


None  of  these  molecules  are  known  experimentally,  yet  the  cyclic  and  polymeric  forms 
conform  to  established  bonding  patterns  with  B  and  C  being  tricoordinate,  with  one  BN 
bond  (1.265  A  being  the  triple  bond  typical  of  iminoboranes  alternating  with  a  longer  single 
BN  bond  (1.482- A),  and  similar  alternation  of  CN  between  single  (1.372  A.)  and  triple 
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bonds  (l.lSrA)  forms.  Besides  being  energetic,  this  polymer  is  likely  to  offer  an  interesting 
semi-conductor  material. 

5.  Other  work  with  metastable  molecules  has  focused  on  tetrahedral  N4  and  the  octahedral 
nitrogen  analog  of  cubane,  [6].  Both  are  found  to  be  metastable;  to  be  symmetry 
(Woodward-Hoffman)  forbidden  from  decomposition  into  N2  molecules;  and  to  be  very 
energetic.  Isp=489  for  N4  and  531  for  Ng.  The  latter  has  been  identified  as  the  “holy  grail” 
of  metastable  molecules,  with  concerted  synthetic  effort  being  direct  toward  it.  We  have 
provided  MBPT(2)  vibrational  spectra  to  facilitate  identification  of  these  species. 

Spectroscopic  Characterization  of  Oxides 

1.  Synthesis  and  detonation  of  energetic  materials  involves  several  critical  oxidation  processes 
and  associated  intermediates.  Carbonyl  oxide. 


1.078 


Figure  7  Carbonyl  Oxide 


and  its  cyclic  isomeric  form,  dioxirane. 


118.8 


Figure  8  Dioxirane 


are  two  key  compounds  in  such  processes.  Carbonyl  oxide  has  never  been  observed 
experimentally,  though  it  is  one  of  the  most  discussed  compounds  awaiting  detection.  We 
performed  high  level  coupled-cluster  CCSD(T)  calculations  to  characterize  carbonyl  oxide 
[10].  These  studies  employed  recently  developed  open-shell  analytical  gradients  methods 
for  CCSD(T),  without  which  structures  and  vibrational  frequencies  would  be  difficult  to 
obtain.  The  AHf°(298)  for  carbonyl  oxide  is  found  to  be  30.2  kcal/mol,  while  the  barrier 
to  isomerization  to  dioxirane  is  19.2  kcal/mol.  CCSD(T)  vibrational  frequencies  of  both 
species  are  presented  to  facilitate  identification  along  with  the  O  isotope  shifts.  Shifts  as 
high  as  45cm~^  permit  experimental  discrimination  between  the  two  forms.  Carbonyl  oxide 
is  found  to  be  far  more  zwitterionic  than  lower  level  theoretical  studies  would  suggest,  in 
line  with  the  viewpoint  of  synthetic  xhemists. 
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2. 


In  the  B+cd2  rdhction,  the  major  product  appeared  to  be  an  unknown  OBCO  species. 
To  identify  it  in  matrix  isolation  requires  the  application  of  predictive  ab  initio  correlated 
methods  to  interpret  the  observed  spectra  (L.  Andrews).  We  applied  coupled-cluster  theory 
(CCSD)  and  MBPT(2)  to  obtain  the  structure  and  vibrational  spectra  for  OCBO  [13].  Five 
potential  structures  were  considered  (Fig.  9). 
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Figure  9  Calculated  structural  isomers  of  BCOj. 


All  were  ruled  out  by  the  computed  vibrational  spectra  except  c  and  e,  with  structure  e  being 
more  stable  by  7  kcal/mol.  The  CCSD  results  for  frequencies  fall  within  5%  of  the  observed 
and  predict  the  isotope  shifts  to  within  an  RMS  error  of  3.4  cm~^  making  assignment  of 
the  OBCO  spectrum  to  the  bent  form  conclusive. 

3.  Similar  applications  have  been  made  for  BO2,  where  there  are  symmetry  breaking  compli¬ 
cations  for  the  theory,  and  are  being  made  for  other  boron  oxides  to  help  resolve  the  lack  of 
experimental  energetic  and  spectroscopic  information  for  the  oxidation  products  of  boron. 


Methodological  Advances 


The  ability  to  provide  accurate  energetic  and  spectroscopic  results  for  candidate  HEDM 
molecules  rests  upon  new  developments  in  ab  initio,  electronic  structure  theory.  Several  such 
advances  have  been  accomplished  in  the  last  period  of  this  grant.  All  are  incorporated  into  the 
ACES  n  program  system  [4]  which  is  now  widely  used  throughout  the  world,  including  industry 
(Lubrizol,  Dupont,  Ford,  etc.)  and  other  government  agencies  (NASA  Ames,  DOE’s  Battelle, 
Pacific  Northwest  Laboratory,  etc.). 

1.  Open  shell  molecules  are  frequently  more  difficult  to  describe  accurately  than  those  where 
all  electrons  are  paired.  The  most  widely  used  procedure  for  open  shells  is  to  start  with 
an  unrestricted  Hartree  Fock  (UHF)  reference  and  add  electron  correlation  via  MBPT  (also 
known  as  MP)  theory..  However,  for  many  problems,  UHF  solutions  suffer  from  spin- 
contamination  as  the  UHF  function  is  not  an  eigenfunction  of  spin.  For  some  problems, 
a  restricted  open  shell  Hartree  Fock  (ROHF)  reference,  which  is  a  spin  eigenfunction,  is 
a  better  choice.  However,  because  the  ROHF  method  uses  different  Fock  operators  in  its 
determination,  one  for  doubly  orbital  orbitals,  and  one  for  singly  occupied  orbitals,  it  was 
not  clear  how  to  use  an  ROHF  reference  in  MBPT  calculations.  We  solved  this  problem  to 
provide  .a  non-iterative  ROHF-MBPT  method  for  second,  third  and  fourth  order,  [2],  that  is 


8 


March  19,  1996 


no  more  expensivfe  than  the  usual  UHF-MBPT  calculations,  but  one  that  offers  vastly  better 
convergence  to  the  exact  result  for  spin  contaminated  cases.  This  is  illustrated  in  Fig.  10, 
^  for  the  CN  radical. 


450 


Level  of  Theory 


Figure  10  Electron  affinity  of  CN  radical.  (O)  ROHF;  (*)  UHF. 


Note  that  at  the  ROHF-MBPT(2)  level,  there  is  good  agreement  with  infinite  order  CCSD 
results,  while  UHF-MBPT(2)  is  far  from  the  converged  answer.  Furthermore,  finite  order 
MBPT  failed  to  rectify  this  failing  of  the  UHF,  spin  contaminated  reference. 

2;  To  locate  critical  points  on  an  energy  surface,- E(R),  as  minima  or  transition  states,  requires 
the  facile  evaluation  of  the  forces  on  the  atoms  in  a  molecule,  VE(R).  Since  molecules  have 
3N-6  degrees  of  freedom  where  N  is  the  number  of  atoms,  repeated  calculation  of  E(R)  in 
all  degrees  of  freedom  is  a  hopeless  procedure  for  most  polyatomic  molecules.  However,  by 
also  evaluating  VE(R)  analytically,  in  about  the  same  time  as  E(R)  itself,  it  is  comparatively 
easy  to  locate  the  minimum  energy  geometry  where  the  atomic  forces  vanish.  We  provided 
the  initial  ROHF-MBPT(2)  analytical  gradients  under  this  grant,  [2,  9]  and  used  it  in  several 
applications  to  transient  molecules  like  FCS. 

3.  The  other  primary  methodological  advances  achieved  two  years  ago  is  the  general  formula¬ 
tion  and  implementation  of  the  equation-of-motion  (EOM-CCSD)  method  for  excited  states, 
UV- visible  spectra,  and  associated  properties.  This  is  a  very  ground  method  that  shows 
promise  of  being  the  method  of  choice  for  excited  states.  The  theory  [11,  12,  16,  23]  and 
several  applications  [14-16,  22]  have  been  recently  reported. 

The  basic  idea  of  EOM-CC  is  that  most  of  the  dynamical  correlation  in  a  ground  (or 
reference)  state  is  quite  similar  to  that  in  an  excited,  ionized,  or  electron  attached  state;  but 
the  non-dynamical  correlation  is  quite  different,  and  any  suitable  method  must  account  for  this 
difference. 
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B.  Equation-of-Motion  Coupled  Cluster  Method  for  Excited  States 

In  the  equation-of-motion  (EOM)-CC,  the  reference  state,  which  is  not  necessarily  the  ground 
state,  is  described  by  CC  theory.  The  fundamental  idea  of  EOM-CC  is  to  describe  the  Ath  excited 
state  wavefunction  as  being  created  from  the  ground  state  by  using  an  excitation  operator  Rk, 
defined  as  a  sum  of  single,  double,  etc.  excitation  operators: 

Rk  =  '^rm{k)q}n  =  ro{k)  +  ...  (1) 

m  i,a  *<•? 

i.e.  a  Cl-like  description.  Amplitudes  are  to  be  determined  via  the  EOM-CC  equations. 
In  this  way,  most  of  the  correlation  effects  that  are  common  to  the  ground  and  excited 
states  are  conveniently  handled,  allowing  the  excited  state  description  to  primarily  focus  on 
the  essential  differential  correlation  effects.  The  Schrodinger  equation  for  the  excited  state 
Hk  =  H-  (0|ff|0)) 


can  then  be  written 

HNRkl^o)  =  HNRke^lO)  =  AEkRke^\0).  (2) 

Multiplying  the  reference  state  solution  i?jv|^o)  =  ^Eo\9o)  by  Rk,  and  subtracting  from 
Eqn.  (2)  yields  the  EOM, 

[iyiv,i?fc]l^o>=‘^it-Rfe|^o)  (3) 

where  Uk  =  AEk  •”  AjE^o  is  the  excitation  energy.  The  operators,  and  Rk  commute. 
Multiplying  Eqn.  (3)  by  e~^  gives 

•  {e-^HNe^)Rk\0)=HNRk\0)  =  ^EkRk\0). 

The  similarity  transformed  Hamiltonian,  Hn,  is  manifestly  connected  and  exactly  terminates  after 
four  T  operators  regardless  of  the  number  of  electrons.  Then  with  AEo  being  the  correlation 
correction  in  the  reference  state,  (^^iv)opg„  =  Hn  ~  AEq,  we  have 

(HN),p,^Rk\0)=<^kRk\0),  (4) 

and  there  can  be  no  contributions  from  the  RkH^  part  of  the  commutator  in  Eqn.  (3),  since  no 
contractions  between  Rk  and  iy^y  are  possible.  In  this  way,  i/jv  is  treated  as  a  generalized,  second- 
quantized  operator  that  can  be  written  as  a  sum  of  one-body,  two-body,  and  higher  interactions: 

Rn  =  fps{ph}  +  Y\  -K  i  +  ...  (5) 

«,f,u 

The  /  and  (||)  integrals  are  modified  or  “dressed”  versions  of  ordinary  Fock  matrix  elements 
whose  specific  expressions  in  terms  of  intermediates  composed  of  two-electron  integrals  and  T 
amplitudes  taken  from  the  reference  state  CC  calculation  are  presented  elsewhere  [11].  For 
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calculations  limited  to  single  and  double  excitations,  i.e.  EOM-CCSD,  nothing  contributes 
beyond  the  three-body  terms. 

From  the  definition  it  is  clear  that  Hn  is  a  non-Hermitian  operator,  so  it  has  left  L  and  R 
eigenvectors,  __ 

{Q\Lk{HN)oj,en='^k{m  (6) 

with  the  same  eigenvalues.  {L*}  and  {Rk}  form  a  set  of  biorthogonal  eigenvectors  {{Lk\Rk)  = 
XkSki  and  we  are  free  to  choose  Ajt  =  1).  Both  Lk  and  Rk  are  required  to  evaluate  transition 
moments  and  other  properties  [11]. 

Eqns.  (4)  and  (6)  can  be  converted  to  a  non-Hermitian  Cl-like  eigenvalue  problem 

HnRh  =  ^kRk, 

LkHN  =  Lk(^k 


Rk  and  Lk  represent  the  coefficients  for  the  chosen  excitations.  Eigenvalues  and  eigenvectors 
of  Hat  can  be  determined  using  a  non-Hermitian  modification  of  Davidson’s  method. 

Since  the  CC  amplitudes  are  solutions,  {Hfi)Qp=Q,  the  first  column  of  (Hjv)  consists  of  zeros, 
which  means  the  excited  state  solutions  do  not  interact  with  Ae  ground  state,  permitting  the 
excitation  energies  to  be  directly  evaluated  from  the  sub-matrix  Notice  ro(K)  is  irrelevant 

for  the  energy  determination.  The  left  eigenvector  for  the  referencee  state  is  Lq  =  1 4-  A,  which 
exactly  comes  from  CC  analytical  gradient  theory,  while  Lk  =  Ajt,  k  >  0.  Consequently, 
EOM-CC  excited  states  can  be  obtained  as  a  byproduct  of  a  CC  analytical  gradient  program 
[11,  12,  16]. 

EOM-CCSD  is  clearly  exact  (i.e.  the  full  Cl)  for  two  electrons,  but  as  will  be  pointed  out 
below,  results  are  much  better  for  states  dominated  by  single  excitations.. 

In  the  last  couple  of  years  under  this  grant  we  have  made  the  initial  general  implementation 
of  the  EOM-CCSD  method  for  electronic  excited  states  and  associated  UV-visible  spectroscopy. 
To  illustrate  the  accuracy  of  the  method  we  can  show  some  comparisons  to  the  “exact”  full  Cl 
results.  The  latter  are  only  possible  for  small  systems  such  as  the  Be  atom  and  the  CH"^  molecules, 
but  they  provide  useful  benchmarks  for  methods  like  EOM-CCSD  that  can  be  applied  to  much 
larger  systems.  Results  for  Be  are  show  in  Tables  I  and  2. 
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Table  1.  Be  singlet  excitation  energies  (in  eV). 
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Table  2.  Be  triplet  exatation  energies  (in  eV) 


2s2p 
2s3s 
2s3p 
2p2  3p* 

2s3d 

2s4s^S 

2s4p^P° 

2s4d^D 

2s5s 

2s5p 

2s6s 


"closed-shell" 

open-shell 

mr-ccsd 

EOM-CCSD* 

E0M-CCSD‘ 

2.80 

2.729 

2.723** 

6.49 

6.447 

6.433** 

7.35 

7.301 

7.290 

7.510 

7.414 

■1 1  ■  « 

7.78 

7.748 

7.730 

8.01 

7.991 

7.973 

8.34 

8.278 

8.267 

8.47 

8.456 

8.437 

8J8 

8J66 

8.548 

8.70 

8.691 

8.681 

8.92 

8.893 

8.874 

deviation 

0.038 

0.013__ 

0.009 

fuua 

2.733^ 

6.444 

7.295 

7.423 

7.741 

7.985 

8.272 

8.449 

8J60 

8.686 

8.886 

0.000 


EXPT 

2.125 

6.457 

7.303 

7.401 

7.694 

7.998 

8.283 

8.423 

8.556 

8.688 

8.823 


0.023 


TViic  Ti/nrV  Used  2s^  ground  state  reference.  .  ^  Sno 

‘  Po, >P s.a«. use. 


reference. 
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The  comparison  between  the  EOM-CC  and  full  Cl  numbers  seen  in  Tables  I  and  2  is 
remarkable,  though  Be  with  only  two  valence  electrons  is  not  the  most  challenging  example. 
With  the  exception  of  the  double  excitation  2s^  2p^,  the  EOM-CCSD  energies  are  all  within 

0.03eV  from  the  full  Cl  values.  Most  excitations  are  within  0.0  leV.  Even  the  double  excitation 
value  is  within  O.leV.  We  will  discuss  improvements  for  states  dominated  by  double  excitations 
below. 

Despite  ground  state  Be  being  a  closed-shell  atom,  to  demonstrate  application  to  open  shells 
we  can  still  calculate  “excitation”  energies  using  an  excited  state  as  the  reference.  Excitation 
energies  using  the  2s2p  and  2s3s  excited  states  as  the  couple-cluster  reference  states 
relative  to  a  UHF  single  reference  (an  ROHE  reference  causes  only  a  slight  difference).  These 
appear  in  Table  2  in  the  column  labeled  “open-shell  EOM-CC.”  These  numbers  compare  very  well 
with  the  ordinary  “closed-shell  EOM-CC”  results.  On  average  the  “open-shell”  numbers  show 
a  slightly  smaller  average  error  (0.009)  that  is  primarily  caused  by  the  superior  description  of 
the  2p^  double  excitation  state,  as  it  is  a  single  excitation  relative  to  the  2s2p  ^P°  reference. 
Clearly,  the  energies  calculated  from  different  reference  states  are  consistent,  providing  great 
flexibility  in  their  choice. 

In  terms  of  the  left  (^|  and  right  |^)  eigenvectors,  the  definition  of  the  dipole  strength  is 

dk. = (7) 

from  which  the  oscillator  strength  is 

fi  = 

Note  the  transition  moment  is  nOt  well  defined,  although  the  actual  observable,  the  dipole 
strength,  is.  In  this  way,  Eqn.  (7)  provides  a  Cl-like  approximation  to  this  quantity. 

As  a  first  test  of  the  accuracy  of  oscillator  strengths  obtained  within  the  framework  presented 
here,  calculations  have  been  carried  out  on  the  two  systems  where  full  configuration  interaction 
(FCi)  results  are  available.  Table  3  documents  excitation  energies  and  dipole  strengths  for  all 
optically  allowed  transitions  that  were  addressed  in  the  FCI  calculations.  Also  included  are 
approximate  excitation  levels  (AEL)  that  act  as  a  measure  of  the  number  of  electrons  that  are 
excited  from  the  CCSD  ground  state  to  form  the  excited  states.  The  AEL  is  defined  by  the 
equation 

AEL^\Y.  (8) 

P 

where  both  the  ground  and  excited  reduced  density  matrices  are  expressed  in  the  quasiparticle 
(approximate  natural  orbital)  basis  that  diagonalizes  the  ground  state  CC  density  [11].  It  should 
be  clear  to  the  reader  that  any  AEL  significantly  larger  than  unity  should  be  regarded  with 
suspicion. 

For  beryllium,  excitation  energies  and  dipole  transition  strengths  obtained  in  the  present 
calculations  are  in  nearly  perfect  agreement  with  the  FCI  results.  Nevertheless,  since  the  EOM- 
CCSD  method  is  exact  for  two  electrons,  and  Be  (Is^  2s^)  has  a  separated  pair  structure,  the 
beryllium  atom  is  an  “easy  “  problem.  The  processes  studied  here  are  essentially  single-electron 
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.  / 

transitions  (all  AEL)  values  are  near  1.05)  and  the  correlation  effects  in  this  system  are  dominated 
by  the  two  valence  electrons. 

A  more  difficult  test  is  provided  by  the  CH"^  molecule,  which  contains  more  electrons  and 
exhibits  strong  nondynamical  electron  effects  in  its  ground  electronic  state  via  interaction  of 
the  l(7^2(r^3(7^  and  la^2cr^l7r^  electronic  configurations.  Due  to  this  important  configuration 
mixing,  one  expects  some  of  the  low-lying  electronic  states  to'have  appreciable  double-excitation 
character,  as  confirmed  by  the  AEL  values  listed  in  the  rightmost  column  of  Table  3.  It  can 
be  seen  that  the  closest  coincidence  of  FCI  and  EOM-CCSD  results  for  excitation  energies  and 
oscillator  strengths  occurs  for  the  transitions  near  3.2  and  13.5  eV,  which  both  have  associated 
AFT,  values  below  1.1.  Differences  between  FCI  and  EOM-CCSD  dipole  strengths  for  transitions 
to  the  third  and  second  state  are  larger  [error  of  5.26%  (AEL=1.13)  and  17.69% 

(AEL= 1.240,  respectively].  The  quantitative  agreement  between  FCI  and  EOM-CCSD  for  the 
dipole  strength  of  the  nearly  pure  double  excitation  at  8.549  eV  is  fortuitous. 

While  the  limited  body  of  FCI  excitation  energies  and  oscillator  strengths  in  the  literature 
prevents  us  from  drawing  many  conclusions  from  the  present  set  of  calculations,  it  appears  that 
the  EOM-CCSD  method  gives  reliable  values  for  dipole  strengths  and  excitations  energies  (and, 
therefore,  oscillator  strengths)  for  electronic  transitions  dominated  by  single  excitations.  More 
recently,  we  have  considered  effects  of  triple  excitations  [20],  which  alleviate  most  of  the  failing 
for  doubly  excited  states.  We  are  currently  using  EOM-CCSD  to  obtain  potential  energy  surfaces 
for  Li  embedded  in  a  cryogenic  H2  lattice  to  help  model  such  embedded  solids. 
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Tible  3.  EOM-CCSD  and  foU  conliguratioB  interaclioo  Mdtatioii  ener^es  and  dipde 
Sangte  f.^  the  loweat  optically  atiowed  ttamitiooa  in  ClT  and  fl.e  berylhnn.  atom.  AUo 
listed  are  the  approximate  exatation  levels  (AEL)  for  each  process. 


ip(2s2p) 

ip(2s3p) 

ip(2s4p) 


5.318 

5.314 

3.564 

3.549 

1.07 

7.465 

7.460 

0.021 

0.022 

1.06 

8306 

8300 

0.000  102 

0.000  121 

1.05 

9.109 

8.549 

0.025 

0.025 

1.96 

13.580 

13.525 

1.073 

1.080 

1.03 

17.315 

17.217 

0.720 

0.684 

1.13 

3.261 

3.230 

0.095 

0.090 

1.03 

14.454 

14.127 

0.692 

0.-588 

1.24 

a.  AoDTOXiinate  excitation  level. 
(1989). 
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Korean  Advanced  Institute  of  Science  and  Technology,  Taejon,  Korea. 

October  1995  —  “An  Introduction  to  Coupled-Cluster  Methods  and  Other  Correlated  Ap¬ 
proaches”,  Postec  University,  Pohang,  Korea. 

October  1995-  “Coupled-Cluster  Theory  as  a  Unified  Approach  to  Molecular  Spectroscopy: 
Vibrational,  Electronic,  and  NMR,”  Korean  Chemical  Society,  Pusan,  Korea 

September  1995  —  “Equation-of-Motion  Coupled-Cluster  Method  for  Excited,  ionized  and 
Electron-Attached  States”,  Molecular  Quantum  Mechanics:  Methods  and  Applications  Confer¬ 
ence,  in  Memory  of  S.F.  Boys  and  in  Honor  of  I.  Shavitt,  University  of  Cambridge,  Cambridge, 
England. 

August  1995  —  “Excited  States  Using  the  Equation-of-Motion  Coupled-Cluster  Method,”  12th 
Canadian  Symposium  on  Theoretical  Chemistry,  University  of  New  Brunswick,  Fredericton, 
Canada 

June  1995  —  “Analogues  of  Tetrahedral  N4  as  Potential  High-Energy  Molecules,”  Air  Force 
High  Energy  Density  Materials  Contractors’  Conference,  Woods  Hole,  MA 

January  1995  —  Air  Force  Office  of  Scientific  Research  HEDM  Theory-Synthesis  Workshop, 
Palmdale,  CA  ■  • 

October  1994  —  “Unknown  Energetic  Molecules  That  Should  Fxist:  A  Qaantum  Chemical 
Investigation,”  46th  Southeastern  Regional  American  Chemical  Society  Meeting,  Birmingham, 
AL 

October  1994  —  “Coupled  Cluster  Theory  as  a  Unified  Approach  to  Molecular  Spectra,”  Work¬ 
shop  on  Non-Perturbative  Many  Body  Methods  (From  Quantum  Field  Theory  to  Chemistry),” 
Bad  Honnef,  Germany 

June  1994  —  “Coupled  Cluster  Theory  as  a  Unified  Approach  to  Molecular  Spectra,”  8th 
International  Congress  of  Quantum  Chemistry,  Prague,  Czech  Republic 

June  1994-  “Properties  in  Coupled  Cluster  Theory,”  8th  International  Congress  of  Quantum 
Chemistry  Satellite  Meeting,  Bratislava,  Slovakia 

June  1994  —  “New  Developments  in  Theory  and  Their  Role  in  Studying  High  Energy  and 
Density  Matter,”  1994  Air  Force  High  Energy  Density  Matter  (HEDM)  Contractors  Meeting, 
Crystal  Bay,  NV 
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May  1994  —  “Excited  States  in  Coupled-Cluster  Theory:  Valence,  Rydberg  and  Core  Excita¬ 
tions,”  2nd  Canadian  Computational  Chemistry  Conference,  Kingston,  Ontario,  Canada 

December  1993  —  Symposium  on  Molecular  Interactions,  Aarhus  University,  Aarhus,  Denmark 

August  1993  —  206th  National  Meeting  of  the  American  Chemical  Society,  Chicago,  Dlinois 

June  1993  —  1993  Air  Force  High  Energy  Density  Materials  Contractors  Conference,  National 
Academy  of  Sciences,  Woods  Hole,  MA. 

June  1993  -  First  Congress  of  the  International  Society  for  Theoretical  Chemical  Physics, 
University  of  Girona,  Girona,  Spain 

May  1993  -  11th  Robert  S.  Mulliken  Lecturer,  University  of  Georgia 

November,  1992  —  Workshop  on  Current  Trends  in  Computational  Chemistry,  Jackson  State 
University,  Jackson,  MS 

July,  1992  —  International  Academy  of  Quantum  Molecular  Sciences,  Menton,  France 

University  CoUoquia  Presented: 

October  1995  —  Lehigh,University,  Bethlehem,  PA 

March  1995  —  Iowa  State  University,  Ames,  lA 
September  1994  —  Ohio  Supercomputer  Center,  Columbus,  OH 
June  1994  —  Edtvos  Lorand  University,  Budapest  Hungary 
-  April  1994  —  Florida  State  University,  Tallahassee,  FL 
February  1994  — Vanderbilt  University,  Nashville,  TN 
July,  1993  —  Freie  Universitat,  Berlin,  Germany 
May,  1993  —  The  Mulliken  Lecturer,  University  of  Georgia,  Athens,  GA 
April,  1993  —  Ohio  State  University,  Columbus,  OH 
March,  1993  —  Vanderbilt  University,  Nashville,  TN 
December,  1992  —  University  of  Virginia,  Charlottesville,  VA 
March,  1992  —  The  Emerson  Lecturer,  Emory  University,  Atlanta,  GA 
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Postdoctoraf  Researchers  Partially  Supported  by  Grant  No. 
AFOSR-F49620-92-J-01 41 

1.  Don  Comeau 

2.  Peter  Szalay 

3.  John  Watts 

4.  Hideo  Sekino 

5.  Jurgen  Gauss 

6.  Marcel  Nooijen 

7.  Anna  Baikova 

8.  Miroslav  Urban 

9.  Stanislaw  Kucharski 

10.  Nevin  Oliphant 

Graduate  Students  Partially  Supported  by  Grant 
No.  AFOSR-F49620-92-J-0141 

1.  David  Bemholdt  —  American  (Ph.D.  1993) 

2.  Sullivan  Beck  —  American 

3.  S.  Ajith  Perera  —  Sri  Lankan  (Ph.D.  1996) 

4.  Lynn  Salemi  —  American 

5.  Seven  Gwaltney  —  American 

6.  Ren6  Mattie  —  American  (Ph.D.  1995) 

Professional  Honors 

Bartlett,  Rodney,  J,  Fellow,  International  Academy  of  Quantum  Molecular  Sciences,  Menton, 
France,  1991 

Bartlett,  Rodney,  J,  18th  most  cited  chemist  in  the  world  for  the  1984-1990  period.  Science 
Watch  Magazine,  May,  1992. 

Bartlett,  Rodney,  J,  American  Physical  Society,  Fellow,  November,  1989 
Bartlett,  Rodney,  J,  Guggenheim  Fellow,  1986-1987 
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Bartlett,  Rodney,  J,  "Chairman,  Subdivision  of  Theoretical  Chemistry,  American  Chemical 
Society,  August,  1987 

Bartlett,  Rodney,  J,  Mulliken  Lecturer,  University  of  Georgia,  1993 
Bartlett,  Rodney,  J,  Cherry  Emerson  Lecturer,  Emory  University,  1992 
Bartlett,  Rodney,  J.  Fulbright  Distinguished  Lecturer,  Univ.  of  Arkansas,  1985. 


